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We develop further the approach of Hubbard and Schofield (Phys.Lett., A40 (1972) 245), which 
maps the fluid Hamiltonian onto a magnetic one. We show that all coefficients of the resulting effec- 
tive Landau-Ginzburg- Wilson (LGW) Hamiltonian may be expressed in terms of the compressibility 
of a reference fluid containing only repulsive interactions, and its density derivatives; we calculate 
the first few coefficients in the case of the hard-core reference fluid. From this LGW-Hamiltonian 
we deduce approximate mean-field relations between critical parameters and test them on data for 
Lennard- Jones, square-well and hard-core- Yukawa fluids. We estimate the Ginzburg criterion for 
these fluids. 



The modern theory of critical phenomena based on 
the Renormalization Group (RG) technique has demon- 
strated an impressive success in a variety of fields [Q. 
However, most of the studies using the RG-approach 
have been addressed to the criticality of the Ising-like 
systems. Some effective computational schemes based 
on the RG-technique have been developed for fluids [H, 
but they are not as convenient for the general analysis 
of criticality as those based on the so-called Ginzburg- 
Landau- Wilson (LGW) Hamiltonian. Effective magnetic 
Hamiltonians for fluids have been derived in a variety 
of ways. In ref. |3j the fluid Hamiltonian was reduced 
to the magnetic one by means of coarse-graining. In 
the coefficients in an effective LGW-Hamiltonian were 
obtained by comparing the critical amplitudes for the or- 
der parameter, compressibility, correlation length, etc., 
calculated within generalized mean-spherical approxima- 
tion, with those derived from the LGW-Hamiltonian. In 
1^ the functional generalization of the Mayer expansion 
for the single-component fluid was used; the attractive 
interactions were treated on the second-virial level, and a 
few different approximation were adopted for the repul- 
sive (hard-core) contribution to the free energy. Within 
these approximations Fisher and Lee evaluated the coef- 
ficients for the effective LGW-Hamiltonian for a single- 
component fluid 1^; somewhat different approximations 
were employed to derive the LGW-coefficients for the re- 
stricred primitive model of electrolytes Q . 

Hubbard and Schofield Q derived the effective LGW- 
Hamiltonians for fluids by an exact mapping, based on 
the transformation of variables. Although they did not 
compute the coefficients of the effective Hamiltonian, 
they argued that fluids belong to the Ising universality 
class. In the present study we develop this approach fur- 
ther, showing in particular that all the coefficients of the 
effective LGW-Hamiltonian may be expressed in terms of 
the known properties of the reference (hard-core) system: 
the compressibility and its density derivatives. We find 
explicit expressions for the first few coefficients. Apply- 
ing the mean-field conditions for the critical point of the 



effective magnetic Hamiltonian we formulate simple rela- 
tions between some critical parameters and check them 
for some fluids; for these fluids we also estimate the pa- 
rameter tq of the Ginzburg criterion. 
We start from the fluid Hamiltonian, 

= E ^("^^^ ) - E ^(^y ) = Hii + HA (1) 

here (f){r) denotes the repulsive part of the interparti- 
cle interaction potential, while —v{r) denotes the at- 
tractive part; {r*,} are coordinates of the particles and 
rij ^ fi — fj. The attractive part of the Hamiltonian, 
Ha, may be written in terms of the Fourier compo- 
nents of the density fluctuations, Uj^ — -^J2f=i 
and the Fourier transform of the attractive potential, 
Vk^j v{r)e-'^^dr, as [|§ 

= E + (2) 

fe 

where N is the number of particles, is the volume of 
the system, and summation over the following set of k 
is implied: fc; = '^ni with L = Q}/'^, I = x,y,z, and 
ni = 0, ±1, ±2, . . .; the thermodynamic limit i — > oo is 
assumed. Let be the chemical potential of the system 
with the total Hamiltonian (|l|), and /i/j the chemical po- 
tential in the reference system having the Hamiltonian 
Hn, which includes only repulsive interaction. Then the 
grand partition function, S, may be expressed in terms 
of that of the reference fluid, Eji, as |^ 

S = Sfl^exp|/?/i'iV-t-i/3^Wfcn£n_,-|^ (3) 

where fi' = fi — /i^ -I- ^v{Q), and {■ ■ ■) fi denotes an av- 
erage over configurations of the reference system, at the 
temperature T with chemical potential /i^. Note that 
the reference system (with only repulsive interactions) 
does not have a liquid-gas transition; i.e. its grand 
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partition function, Sjj, is regular in the vicinity of the 
critical point of the fluid of interest. Following Hub- 
bard and Schofield we use the identity exp(ia^a;^) = 
(27ra^)^^/^ exp{—^y'^/a'^+xy)dy, and after some al- 
gebra obtain the ratio Q = S/S/j of the grand partition 
functions [||: 
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here integration is to be performed under the restric- 
tion cr_^ = at (the complex conjugate of cr^); ctq = crfc 
for fc = 0, and we omit a factor regular at the criti- 
cal point since it does not affect the subsequent analy- 
sis. Applying the cumulant theorem Q to the factor 

exp l^jT (Tj:n_^|^ one obtains ||^: 
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n > 3 



cR. 



Here p = iV/fi is the fluid density and {■ ■ ■)cr denotes the 
cumulant average calculated in the reference system. In 
Eqs.(||) Q has been written in the same way as the parti- 
tion function for a magnetic system having an Ising-like 
Hamiltonian: are the Fourier components of the "spin 
fleld", cr(r), and h' is the "magnetic field". 

The coefficients of the effective Hamiltonian depend 
on the correlation functions of the reference fluid having 
only repulsive interactions. Using the deflnitions of the 
particle correlation functions of fluids and definitions 
of the cumulant averages j^, one can directly evaluate 

(nj:^ . . .nj: ^ , and thus the coefficients u^(fci, . . . , fc„). 

It is straightforward to show that may be expressed in 
terms of the Fourier transforms of the correlation func- 
tions hi, /i2 . . ., hn of the reference system, defined as 
hi{ri) = S{fi), h2iri,f2) = .92(^1, ^^2) - 1, ^13(^1, r2, fa) = 
53(^1,^2,^3) -32(^1,^2) -52(^1,^3) -52(^2,^*3) + 2, etc. 
where gi{fi, . . . ,ri) are /-particle correlation functions 
In particular, the first few coefficients read: 
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- {l + ph2{ki) 
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(^1,2,3 P i 1 + P h2{ki) + h2{k2) + ^12(^3) 



+p2^3(fci;fc2)} 



(7) 



where hi are the Fourier transforms of hi , and we use the 
shorthand notation: (5i,2,...,r! = ^j^+g^ +k o/"''' 

Now we consider the small-fc expansion of the co- 
efficicients. First we note that the function h2{k) 
may be expressed in terms of the Fourier transform 
of the direct correlation function, C2(fc), as h2{k) = 

C2{k) / ^1 — pc2{k)j , and its zero-fc value in terms of the 
isothermal compressibility, XR = P^^ {(^p/ (^Pr) (where 
Pr is the pressure of the reference fluid) as 1 -I- ph2 (0) = 
pksTxR. = zq- Using the expansions — v,^— Wq fc^ + • • ■ 

and C2(fc) = £2(0) - C2(0)"fc2 H , one obtains for 

M2, omitting terms of C'(fc'*), 



u'2 — 5i^2 [a'2 + b^k^ + ■ ■ ■] 
a'2 = {PvaV^ - pzo 



(8) 
(9) 



The LGW-Hamiltonian does not have terms with powers 
of k higher than fc^; moreover, the only term of order k"^ 
reads oc k^a^a_^. Thus, only zero-order terms should be 
kept in the expansion of u'^ for n > 2. Hence we may 
write the contribution of such terms using this approxi- 
mation: 



-(51,2,3 P 



1 + '3ph2i6) + p^hiO) 



(10) 



< - -'^1,2,3,4 P [1 + 7p/i2(0) + 6p^h{0) + P^hiiO)\ , 

etc., where hi{0) = hi{0, ... ,0). There exists a relation 
between successive correlation functions H , 
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from which follows a relation between the functions hi{Q) 
d 



xp^^P%{6) = /3p' I hi{6) + phi+i{6) 
dp 



(12) 



expressing each /i/+i(0) in terms of /i;(0) and its den- 
sity derivative. Using Eq.(|l^ itcratively one flnally finds 
each hi{0) expressed in terms of the reference system 
compressibility XR ^-nd its density derivatives. Explic- 
itly we obtain: 

"3 = '^1,2,3 PcU3 = -(5i,2,3 pza [zo + zi) (13) 

u'4 = (5l,2,3,4 PcW4 = 

= -(^1,2,3,4 PZQ [zf + Zq (zq + 4zi + Z2)] , (14) 

defining m„, where pc is the critical density, zq = pxr/ P 
as before, zi = pdzo/dp and Z2 = p^d'^z^jdf? . In this 
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way one can evaluate all the coefficients of the effective 
LGW Hamiltonian of the fluid and express them in terms 
of the compressibility of the reference system and its den- 
sity derivatives. This solves the problem of finding the 
effective LGW Hamiltonian, provided that the compress- 
ibility of the reference system is sufficiently well known. 

For the reference system with only repulsive interac- 
tions one can often usefully adopt the hard-sphere system 
with an appropriately chosen hard-core diameter ||ll|,|l2| . 
For soft (not impulsive) repulsive forces a simple relation 
^d = dBH = /o" [1 - exp {-(l>{r)/kT)] gives the effec- 
tive diameter d of the hard sphere system correspond- 
ing to a repulsive potential 0(r) that vanishes at r > a. 
For the hard-sphere system one has the fairly accurate 
Carnahan-Starling (CS) equation of state for which 



zo = (1 - 77) (1 + 477-1- V - 4?7^ + V 



4^-l 



(15) 



where rj = ^cPp. The value of £2(0)" may be found from 
the Wertheim and Thiele solution ||l^ for the the direct 
correlation function, which gives 

C2(0)" = - (7rdVl20) (16 - Ih] + At^^) (1 - 77)"^^ . (16) 

To recast the effective Hamiltonian in the conventional 
form, we perform a transformation from the variables 
aj: to "field" variable cr^r). Under this transformation 
integration over the set {crj} becomes "field" integra- 
tion over a{f), and the term ^ k^a^a_^ transforms into 

- {V(j{r)y. Using Pc as a scaling factor for the length 
we finally arrive at the effective LGW-Hamiltonian: 



n 



df 



-ha + —a + —a + —<7 



3! 



4! 



|(v.)^; 

(17) 



where h = {p - ^jlr + v{Q) /2) {v^pc) ^ -I- {p/Pc), 02 = 
{kBTJpcVo)-zo-{p/pc), and M3, U4 are given by Eqs.(|l|) 
and ([l4|). The coefficient 62 reads: 
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X2 



B 



(18) 



where B = 477^(1 - 77)^^(16 - II77 -I- 4772)7(1 + 'i-q -\- Arf' - 
Arf + rfY and constants Ce// and Ag// characterize zero- 
order, Wo, and second-order, Vq = J r'^v{r)dr, moments 
of the attractive potential: 



47rd3 



// = / vir)df^vo 
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(19) 
(20) 



This effective LGW-Hamiltonian may be used for anal- 
ysis of the critical behavior, using for example the RG- 
technique. 



Now we perform a mean-field level (MF) analysis based 
on to the effective Hamiltonian. To do so we first remove 
the cubic term by making the shift, <t ^ <t + W, with a 
chosen to make the cubic term vanish; this leads to new 
coefficients: h — h + a2U3/u4 — u^/3ul, 02 = 02 — 'u§/2it4 
with U3 = 0, U4 = U4 and 62 = 52- Then the MF- 
condition at the critical point, 02 = and h = 0, gives 
the approximate relations 



PcVo 



VoPc 
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6ul 



2U4 



(21) 
(22) 



c.p. 



Eq.(pTl) relates the critical density and the critical tem- 
perature of the system; Eq. (p^) relates the difference of 
the chemical potentials, p — pu at the critical point to 
the critical density. We have tested these relations us- 
ing simulation data for Lcnnard- Jones (LJ) [Q, hard- 
core-Yukawa (HCY) |l|l and square-well (SW) (l|-|l|] 
fiuids. The WGA partition |0] of the LJ-potential, 
ULjir) = 'A^LJ [(o'/7')^2 _ [a/rY\, gives for the attrac- 
tive part. 



v{r) 



-UL.j{r) 



r < 25(7 
r > 25(7 



(23) 



which is perfectly smooth in the core region. This par- 
tition gives the best estimates for the thermodynamic 
functions in the WCA-perturbation scheme The 
repulsive part, (f>{r) = ULj(r) + w(r), is then used to 
find the effective hard-core diameter using the expres- 
sion quoted above. Similarly we use the WCA-partition 
for the SW and HCY fluids. The square- well fluid has 
the interaction potential, (/)(r) — v(r) = -l-oo if r < d, — e 
if (J < r < Ad and if r > Ad; the reference system is 
the hard sphere system with the diameter d. We take the 
attractive part of the potential as v{r) — for r > Ad, 
and v{r) — e for < r < Ad. For the HCY potential, 
4>{r) — v{r) = +00 if r < d, and — ey exp [—K{r — d)] /r 
for r > d; the WCA-partition gives v{r) = ey/d for r < d 
and v(r) = ey exp [—K{r — d)] /r for r > d. 

Table I gives the ratio of the right over the left-hand 
side of the MF Eq.(|2l|) as Wc, and that of Eq.(|2|) as Yc. 
As one can see from the Table, these MF-relations hold 
rather satisfactorily except for HCY-fluids with short- 
ranged attractive potentials, where the MF-description 
and MC results apper to differ. In Table I we also give 
some coefficients of the LGW-Hamiltonian at the critical 
point. 

Using coefficients of the effective LGW-Hamiltonian 
obtained, one can also estimate the Ginzburg parameter 
tq |jl|. This defines the domain of validity of the clas- 
sical critical behavior: the classical description fails for 
r| = \T/Tc - 1| < tg- Following ||^, we write for this 
parameter: 



3 



32 TT- 



(24) 



where 0,2 = a2T. From Eqs.(|l8[) and (24) it follows that 
for the infinitely-ranged Kac-Baker potencial, Xeff 
00, (with eg// oc vq finite, see Eqs.(p^ and ( pO| ) ), 
&2 ^ 00, and thus tg ^ as expected (of. t o 

To estimate a2 we use the MF-condition (pl[) and ap- 

proximate zg • [p/ Pc) + 2ir value at the critical 



point, 



PcVO 



This yields: 



a2 



a2r 



(25) 



kBTc/8ric£ef f ■ Using this value of a2, 
^h, and (Eq) (the coefficients in 



Thus, QfT 

and U4, 62 given by Eqs 
this relation computed at the critical point), we calculate 
an approximate Ginzburg parameter tq for some of the 
SW, HCY and LJ-fluids. The results are given in Table 

As we see in Table I, the derived value of tq is of the 
order of « 10~^ for the most of fiuids studied in com- 
puter simulations and lies within the range of tq values 
predicted in Ref. On the other hand, for the SW- 
fluid with the most long-ranged attraction {Xeff = 3), 
tg ~ 10~^, is much smaller, and may explain the MF- 
like behavior observed iQ. For the LJ-fluid, with an 
attractive potential, oc r~^, similar to that of real fluids, 
tq is of the order of 10~^; real simple fluids are supposed 
to have similar values of tq |20|| . However, it may be 
noted that in MC simulations of the LJ-fiuid (where 
a cut-off of the LJ-potential and tail corrections were 
used) Ising-like behavior was observed to values of |t|, at 
least few times larger than the tq given in Table I. 
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Table I Gives the Ginzgurg criterion, tq, the ratio 
of the right over the left-hand side of Eq.(|2^) as Wc, 
and that of Eq.(|2^) as Yc, and coefficients of the effec- 
tive LGW-Hamiltonian, &2,c, U4,c, 0^2- Here Ag// is the 
effective range of the attractive potential (Ae// = A for 
the SW-fiuid), T* = keT/teff and pi = pcdP are re- 
spectively, the reduced critical temperature and density. 
The critical data for SW-fluid is from [^ , from jl^] (for 
= 1.5), from jl| (for Ae// = 3.0); for the HCY-fluid 
from ITsIl ; for the LJ-fluid from ITll . 
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